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ABSTRACT 

We present an analysis of the broadband radio spectrum (from 22 to 3900 MHz) of the Galactic 
supernova remnant (SNR) HB3 (G132.7+1.3). Published observations have revealed that a curvature 
is present in the radio spectrum of this SNR, indicating that a single synchrotron component appears is 
insufficient to adequately fit the spectrum. We present here a fit to this spectrum using a combination 
of a synchrotron component and a thermal bremsstrahlung component. We discuss properties of this 
latter component and estimate the ambient density implied by the presence of this component to be n 
~ 10 cm -3 . We have also analyzed extracted X-ray spectra from archived ASCA GIS observations of 
different regions of HB3 to obtain independent estimates of the density of the surrounding interstellar 
medium (ISM). From this analysis, we have derived electron densities of 0.1-0.4 f^ 1 ^ 2 cm" 3 for the 
ISM for the three different regions of the SNR, where / is the volume filling factor. By comparing 
these density estimates with the estimate derived from the thermal bremsstrahlung component, we 
argue that the radio thermal bremsstrahlung emission is emitted from a thin shell enclosing HB3. 
The presence of this thermal bremsstrahlung component in the radio spectrum of HB3 suggests that 
this SNR is in fact interacting with an adjacent molecular cloud associated with the HII region W3. 
By extension, we argue that the presence of thermal emission at radio wavelengths may be a useful 
tool for identifying interactions between SNRs and molecular clouds, and for estimating the ambient 
density near SNRs using radio continuum data. 

Subject headings: radiation mechanisms: thermal — radio continuum: general — supernova remnants 
- ISM: individual: HB3(G132. 7+1.3) 



1. INTRODUCTION 

It has been firmly established that the dominant 
emission mechanism from Galactic supernova remnants 
(SNRs) at radio frequencies is synchrotron emission. 
This conclusion has been reached based on the measured 
spectral index a of the observed radio emission (S v ~ 
v "~") from these sources ranges from w 0.3 - 0.7 (| Green! 
2006) and - in some cases - the detection of polarized ra- 
dio emission. Typically, the radio emission from Galactic 
SNRs can be modeled by a single power law: however, in 
the cases of some SNRs, observations over a very broad 
range of radio frequencies reveal a curvature in the spec- 
tra of these sources. Two scenarios have been presented 
in the literature to explain this observed curvature: the 
first scenario invokes spectral index variations due to dif- 
ferent populations of synch rotron emitting electrons as- 
sociated with the SNR (e.g-. lTian fe Leahvl (|2005t )). while 
the second scenario argues for the presence of a thermal 
bremsstrahlung component in the radio emission from 
the SNR in addition to the synchrotron component. This 
thermal bremsstrahlung component is expected to be 
most promi nent for evolved SNRs int eracting with molec- 
ular clouds (|Urosevic et alj |2003a,b: Urosevic fc Pannutil 
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To critically evaluate these two proposed scenarios, 
we present an analysis of the broadband (38-3900 MHz) 
radio spectrum of one Galactic SNR known to feature 
a curvature in its radio spectrum, HB3 (G132.7+1.3). 
This source has been the subject o f extensive radio 
conti n uum, HI and OH obser va tions (|Brown fc Hazard! 
1 19531: [Landecker et all Il987t iRoutledge et alj 11991k 
iKoralesky et all 119981 : iTian fc Leahvl l2005h which have 
revealed a shell-like radio morphology for this source 
with an angular diameter of approximately 80 arcmin- 
utes. The spectral index of this source is a = 0.4 and 
the measured integrated flux density at 1 GHz is 45 Jy 
(|Greenll2006f ): radio observations of this SNR are com- 
plicated by confusing thermal emission from the adja- 
cent HII region W3. Thermal emission from an X-ray- 
emitting plasma located in the inte rior of HB3 was dis - 
covered by Einstein and described bv lLeahv et al.l (|1985h . 
This X-ray emissi on is seen to lie ent irely within the ra- 
dio shell of HB3 (jLeahv et al.l IT985) : this combination 
of a radio shell morphology with a center-filled thermal 
X-ray morphology has led t o the classification o f HB3 as 
a mixed-morphology SNR (|Rho fc Petrel 1 1998h . In this 
paper, we separately model the broadband radio spec- 
trum of HB3 with a thermal bremsstrahlung component 
combined with a synchrotron component: from this mod- 
elling wc obtain an estimate of the density of the ambient 
interstellar medium (ISM) surrounding HB3. We com- 
pare this estimate to ambient density estimates obtained 
from analyzing extracted X-ray spectra of this SNR: we 
then discuss a possible interaction between HB3 and W3. 
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2. MODELING THERMAL BREMSSTRAHLUNG EMISSION 
AT RADIO FREQUENCIES FROM HB3 

As descr ibed in previous works ( Urosevic et al.1 

l2003alTbl ; lUrosevic fc Pannutil I2005T ). thermal 
bremsstrahlung emission may become a significant 
emission process at radio frequencies for evolved SNRs 
which are expanding into regions of the interstellar 
medium (ISM) with enhanced densities such as a 
molecular cloud. The volume emissivity e v (in CGS 
units) of thermal bremsstrahlung emission for an ionized 
gas cloud is proportional to the square of the electron 
(or ion) volume density n, i.e., 

£„[crgs s- 1 cm- 3 Hz" 1 ] = 7 x l(T 38 n 2 T- 1/2 (1) 

(Rohlfs & Wilson 1996), where T is the thermodynamic 
temperature of the ISM. The presence of a significant 
amount of thermal bremsstrahlung emission will pro- 
duce a curvature in the observed radio spectrum of the 
SNR, particularly fo r frequencies of 1 GHz and greater. 
iTian fc Le ahy (2005) showed that spectral flattening ex- 
ists in the radio spectrum of HB3 at higher frequen- 
cies and argued that this flattening originates from syn- 
chrotron emission from different populations of electrons 
with a range of energy (spectral) indices in low, middle 
and high radio frequency domains. For comparison, we 
have modelled the radio spectrum of HB3 using a sin- 
gle synchrotron component combined with an additional 
thermal bremsstrahlung component. For this artificially 
generated spectrum, we assume that the amount of ther- 
mal bremsstrahlung emission is 40% and the amount of 
synchrotron emission is 60% of total emission at 1GHz. 
In addition, we generated an artificial spectrum using 
an approximately "pure" synchrotron component with a 
spectral index a — 0.65. This value for a was obtained by 
fitting the radio flux density from HB3 over the low fre- 
quency range 22 to 178 MHz (we expect that the thermal 
component should account for 10-25% of the observed 
radio flux density over this frequency range) and a ther- 
mal component with a spectral index a — 0.1. Based 
on these parameters, we calculated the frequency spec- 
trum over 10 7 - 10 11 Hz. In Figure 1 (left) we present 
both the resultant spectrum (as plotted with a thick line) 
and the measured flux densit ies for HB3 over thi s fre- 
quency range as compiled by ITian fc Leahvl (|2005D and 
references therein. From inspection of this Figure it is 
clear that the inclusion of a thermal component produces 
a noticeable bend in the radio spectrum which fits the 
observed radio spectrum well. We therefore argue that 
the total radio spectrum of HB3 is successfully modelled 
as the sum of synchrotron and thermal bremsstrahlung 
components. 

Based on the properties of the thermal component used 
to fit the radio spectrum of HB3, we can estimate the am- 
bient density of the ISM into which the SNR is expanding 
as follows. Using values obtained from radio observations 
(such as flux density S„, diameter D and thickness of 
SNR shell s) and assu ming a distance to HB3 of 2 kpc 
(jRoutledge et al.lll991h and an electron temperature T 
= 10 4 K, we use the relation given in Eq uation 1 to cal- 
culate the volume emissivity. Specifically, ITian fc Leahvl 
(2005!) reported Si G Hz = 50 Jy, D — 70 pc. Using s 
= 0.05D as a reasonable estimate for shell thickness of 
evolved SNRs, we obtain e v = 1.67xl0 -37 ergs sec -1 



cm -3 Hz^ 1 . From Equation [T] and assuming that the 
thermal component produces 40% of the total radio emis- 
sion from HB3 at 1 GHz, we estimate that the density of 
the ISM into which HB3 is expanding is w 10 cm -3 . 



3. X-RAY EMISSION FROM HB3 

Because the X-ray emission from HB3 extends over 
such a large angular extent on the sky, pointed obser- 
vations were conducted with the Advanced Satellite for 
Cosmology and Astrophysics (ASCA) of three different 
regions (hereafter referred to as the northern region, the 
central region and the southern region) to ensure that 
virtually all of the X-ray emission from HB3 was sam- 
pled. The bulk of this emission originates from the cen- 
tral region where a prominent ring (appr oximately 35 
arcminutes in diameter) of emission is seen ([Leahy et al.l 
1985). The significant variation in the X-ray brightness 
of HB3 from one portion of the SNR to another clearly 
indicates that a wide range exists in the density of the 
X-ray emitting gas. To help quantify the range of this 
variation, we analyzed extracted spectra of the emission 
from each of the three sampled regions as sampled by 
the two Gas Imaging Spectrometers (GIS) - denoted as 
GIS2 and GIS3 - which were aboard ASCA and collected 
data during these observations. Details of the GIS ob- 
servations of HB3 are provided in Table [T] 

The data reduction process for the two observations 
were conducted using the "XSELECT" program (Ver- 
sion 2.2), which is available from the High Energy As- 
trophysics Science Archive Research Center (HEASARC 
6 ). We extracted both GIS2 and GIS3 spectra using the 
following extraction regions: for the northern region, we 
used a circular region with a radius of ~10.3' (centered 
at RA (J2000.0) 02^ 20 m 06 s , Dec (J2000.0) +63° 17' 
44"); for the central region, we used an elliptical region 
with radii 14.7' x 14.2' (centered at RA (J2000.0) 02 h 
17 m 17 s , Dec (J2000.0) +62° 41' 49") and finally for the 
southern region we used a circular region with a radius 
of ~12.3' (centered at RA (J2000.0) 02 h 18 m 15.3 s , Dec 
(J2000.0) +62° 10' 05"). The spectra were extracted in 
the format of pulse-height amplitude (PHA) files as well 
as images and event files: high bit-rate and medium bit- 
rate telemetry data were used and the standard REV2 
screening criteria were applied. For the analysis of the 
GIS2 and GIS3 spectra, we used the standard response 
matrix files (RMFs) and the corresponding ancillary re- 
sponse files (ARFs) were prepared using the program 
"ASCAARF" ( Version 3.10) fro m the FTOOLS soft- 
ware package 7 (|Blackburnl fl995). Background spectra 
and images for both GIS2 and GIS3 were prepared us- 
ing the FTOOL "MKGISBGD" (Version 1.6). In Figure 
1 (right), we present the total band (0.7-10.0 keV) mo- 
saicked ASCA GIS2+GIS3 image of X-ray emission from 
HB3 with radio (1420 MHz) contours overlaid. 

6 HEASARC is a service of the Laboratory for High Energy 
Astrophysics (LHEA) at the National Aeronautics and Space Ad- 
ministration Goddard Space Flight Center (NASA/GSFC) and 
the High Energy Astrophysics Division of the Smithsonian Astro- 
physical Observatory (SAO). For more information on HEASARC, 
please see http://hcasarc.gsfc.nasa.gov 

7 "ASCAARF" and all of the other programs mentioned 
in this Section are part of the FTOOLS software package. 
For more informatio n on this so ftware package, please see 
http://heasarc.gsfc.nasa.gov/ftools/ 
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To obtain an independent estimate of the ambient den- 
sity of the ISM surrounding HB3, we performed spectral 
fitting on the extracted ASCA GIS spectra for the three 
different regions of this SNR. For each of the three re- 
gions, we simultaneously fit the extracted GIS2 and GIS3 
spectra with a collisional plasma model com ponent with 
solar elemental abundances known as APEC (jSmith et all 
2001) 8 combined with the PHABS component to model 
interstellar extinction. Using this combination of compo- 
nents, we were able to derive statistically acceptable fits 
(x 2 „ ~ 1.2 or less) for the extracted spectra of the north- 
ern and southern regions, with column densities ranging 
from Nh ~ 0.2 - 0.6 x 10 22 cm~ 2 and temperatures 
kT « 0.4 - 0.6 keV. A statistically acceptable fit to the 
extracted GIS spectra of the central region could not be 
obtained with this combination of components. Promi- 
nent emission lines associated with magnesium and sil- 
icon were seen in the extracted spectra, so a collisional 
plasma model component with variable elemental abun- 
dances - VAPEC - was implemented with the abun- 
dances of those two elements allowed to vary. A Gaussian 
component was also added to account for an emission line 
of unknown origin seen near 1.2 keV. A statistically ac- 
ceptable fit still was not obtained, so additional attempts 
to fit the spectra using a second component - either an 
additional APEC component with a higher temperature 
or a power law component - were made. A statistically 
acceptable fit was at last obtained after the addition of 
either of these later two components. The column densi- 
ties and temperatures of the lower temperature thermal 
components are comparable in the two fits (Njj ~ 0.4- 
0.5xl0 22 cm~ 2 and kT ~ 0.2 keV, respectively). The 
temperature of the higher temperature component is kT 
w 2.2 keV while the photon index of the power law com - 
ponent is T w 2.6. We note that iLeahv et ail (|1985f ) 
has already discussed the presence of two thermal com- 
ponents in the X-ray emission from the central region 
of HB3 as revealed by Einstein observations. We also 
note that a gradient is seen in the column densities of 
the three different regions, decreasing in proceeding from 
south to north: this result is consistent with the presence 
of denser material seen in projection toward the southern 
half of HB3. We summarize the results of these spec- 
tral fits in Table O A separate analysis of the ASCA 
data for HB3 has also b een analyzed and discussed by 
lLazendic fc Slanel (pOOl . 

Using the value for the emission measure of the soft 
thermal component as derived in fits to the extracted 
GIS spectra for all three regions, we calculated the 
corresponding electron density of the ambient ISM for 
each region. We set each emission measure equal to 
(10~ 14 /47roP)xn e xnnxfxV, where d is the distance 
to HB3 (assumed to be 2 kpc), n e and nn are the elec- 
tron and hydrogen densities (with n# = 1.2 n e ), / is the 
volume filling factor and finally V is the volume (assum- 
ing that the depth of the emitting region is equal to its 
radius). From this relation, we calculate electron densi- 
ties n e w 0.4 J -1 / 2 cm~ 3 for the central region and n e 
w 0.1 f~ x l 2 cm" 3 for the northern and southern regions 
(see Table [2]). 

4. DISCUSSION 



For comparison to Tian and Leahy (2005): Table 2 
there summarizes the spectral variations: the spectral 
indices range from to 0.9 with most values between 
0.3 and 0.7. The boxes nearest W3 are 9 and 13: they 
both contain flatter indices than 0.3, especially box 13. 
The variable spectral index was suggested to be due to 
multiple populations of electrons since the normal case 
of a single population with steeper index at higher en- 
ergies does not work to explain the radio spectrum (e.g. 
their Fig. 3). We are suggesting here that a more natu- 
ral explanation is synchrotron emission which dominates 
at lower frequencies and bremsstrahlung emission which 
dominates at higher frequencies. This is supported by 
the high thermal electron densities derived from the X- 
ray observations. One can still produce spatial variation 
of observed spectral index caused by the spatial variation 
of relativistic electrons and magnetic field (for the syn- 
chrotron component) and of thermal electrons (for the 
bremsstrahlung component). Here we present only a sin- 
gle fit to illustrate that the general idea of synchrotron 
plus bremsstrahlung emission is viable. 

By comparing the densities implied by our fit 
to the broadband radio spectrum with a thermal 
bremsstrahlung component with the densities implied by 
our X-ray observations, we argue that the density profile 
is such that ambient material is gathered in a thin shell 
on the outer edge of HB3: in fact, the observed thermal 
bremsstrahlung emission at radio frequencies is produced 
from this shell. We also claim that the distribution of this 
material in a shell indicates that HB3 is in fact interact- 
ing with the adjacent molecular cloud W3. There has 
been some debate pre viously in the lit e ratur e regarding 
such an interaction: iRoutledge et al.1 (|1991l ) described 
the detection (through CO observations) of the molecu- 
lar bar feature that is partially surrounded by continuum 
emission from HB3, suggesting that HB3 is indeed inter- 
actin g with this cloud (and , by extension, W3). How- 
ever, |Kor^esk^iL a l- (1998) failed to find shock-excited 
maser emission at the nominal boundary between the 
HII region and the SNR that could clearly be associated 
with the shock of HB3. Our results supports interaction 
between HB3 and W3; we also note that the gradient 
in the column densities inferred by the X-ray observa- 
tions is consistent with the geometry of the known cloud 
and further supports such an interaction. In fact, the 
detection of thermal emission at radio frequencies from 
Galactic SNRs may be a crucial new tool in determin- 
ing whether these sources are interacting with adjacent 
molecular clouds and for estimating the ambient density 
near SNRs using radio continuum data. 
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of the Project No. 146012 supported by the Ministry of 
Science and Environmental Protection of Serbia. 
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Fig , 1. — (left) The broadband radio spectrum of HB3 based on measured flux densities for this SNR as compiled by ITian fc Leahvl 
(2005) as fit with the two-component model (synchrotron and thermal bremsstrahlung). (right) ASCA GIS mosaicked image of HB3: the 
image has been background-subtracted, exposure-corrected and smoothed with a Gaussian of 1 arcminute. The pixel range of the X-ray 
emission is to 1.43X10 -4 counts sec 1 pixel 1 . Radio emission (408 MHz) has been overlaid in white contours: the contour levels are 
60, 65, 70, 75, 80, 85, 90, 95, 100 110 and 120 mjy/beam. The regions of spectral extraction are marked with circles and ellipses. 



TABLE 1 

Summary of ASCA GIS Observations of HB3 (G132.7+1.3) 













GIS2 


GIS2 


GIS3 


GIS3 












Effective 


Count 


Effective 


Count 






Sampled 


Right 




Exposure 


Rate 


Time 


Rate 


Sequence 


Observation 


Region 


Ascencion 


Declination 


Time 


(cts 


Exposure 


(cts 


Number 


Date 


of HB3 


(J2000.0) 


(J2000.0) 


(sec) 


sec^ 1 ) 


(sec) 


see ) 


54009000 


27 August 1996 


South 


02 18 43.7 


+62 08 29 


13406 


0.2591 


13404 


0.2660 


54009010 


28 August 1996 


Center 


02 17 07.5 


+62 42 55 


29720 


0.3806 


31633 


0.4075 


54009020 


29 August 1996 


North 


02 19 53.9 


+63 16 53 


9261 


0.2399 


9260 


0.2585 



Note. — The units of Right Ascencion are hours, minutes and seconds and the units of Declination are 
degrees, arcminutes and arcseconds. Count rates are for the energy range 0.7—10.0 keV. 



Also see Jhttp: / /hea- www. haravrd.edu/APEC 



REFERENCES 



Blackburn, J. K., ASP Conference Series, Vol. 77, 1995, R. A. Shaw, 

H. E. Payne and J. J. E. Hayes, eds., p. 367 
Brown, R.H. & Hazard, C. 1953, MNRAS, 113, 123 
Green, D. A. 2006, A Catalogue of Galactic SNRS (April Version), 

available at http://www.mrao.cam.ac.uk/surveys/snrs/ 
Koralesky, B. et al. 1998, AJ, 116, 1323 
Landecker, T. L. et al. 1987, AJ, 94, 111 
Lazendic, J. &: Slane 2006, ApJ 
Leahy, D. A. et al. 1985, ApJ, 294, 183 
Rho, J. & Petre, R. 1998, ApJ, 503, L167 



Rohlfs, K. & Wilson, T.L. 1996, Tools of Radio Astronomy, 
Springer 

Routledge, D. et al. 1991, A&A, 247, 529 
Smith, R. K. et al. 2001, ApJ, 556, L91 
Tian, W.W. & Leahy, D. 2005, A&A, 436, 187 
Urosevic, D. et al. 2003a, Serb. Astron. J., No. 166, 61 
Urosevic, D. et al. 2003b, Serb. Astron. J., No. 166, 67 
Urosevic, D. & Pannuti, T.G. 2005, Astropart. Phys., 23, 577 



Thermal Radio Emission from HB3? 



5 



TABLE 2 

Summary of Spectral Fits to ASCA GIS Spectra of the Different Regions of HB3 a 







PHABSx(VAPEC 


PHABSx(VAPEC 




Parameter 


PHABSxAPEC 


+APEC+Gaussian) 


+Power Law+Gaussian) 


PHABSxAPEC 


Region 


South 


Center 


Center 


North 


2/2 /T\/~\T^\ 

Xt (X /DOF) 


1.15 (451.33/394) 


1.19 (458.26/386) 


1.21 (467.82/386) 


1.22 (481.18/394) 


Njj (10 cm z ) 


°- 55 -0.15 


°- 44 -0.12 


°- 47 -0.27 


n _. „+o.08 
°- 16 l0.06 


kT low (keV) 


U.J9_ Q 09 


°- 22 -0.04 


n oi +0.11 
°- 21 -0.07 


o K/l+0.05 
°' 54 -0.10 


Mg 


1.0 (frozen) 


2.41^ 


2.9(>1.0) 


1.0 (frozen) 


Si 


1.0 (frozen) 


4-2+1° 


4.6(>2.0) 


1.0 (frozen) 


EM io „ b /(47rd 2 /10- 14 ) (cm- 5 ) 


2xl0- 2 


0.17 


0.19 


3xl0- 3 


n e (cm' 3 ) 


0.14/" 1 / 2 


0.35/" 1 / 2 


0.38/" 1 / 2 


0.08/- 1 / 2 


Egauss (keV) 




1.26 


1.22 




Normalization 




4.1xl0- 4 


2.0X10- 4 




kT high (keV) 




9 9+3.8 
z - z -0.2 






EM W9h b /(47rd 2 /10- 14 ) (cm- 5 ) 




1.3xl0- 3 






r 






2-6+1:! 




Normalization d 






7.4xl0- 4 




Absorbed Flux e (ergs cm -2 sec -1 ) 


4.2xl0- 12 


1.6X10- 11 


1.6X10- 11 


3.5xl0- 12 


Unabsorbcd Flux° (ergs cm -2 sec— 1 ) 


2.4X10- 11 


6.8X10- 11 


7.0X10- 11 


5.9xl0- 12 


Luminosity (ergs sec- 1 ) 


l.lxlO 34 


3.3 xlO 34 


3.4xl0 34 


2.8xl0 33 



a All quoted errors are 90% confidence intervals. 13 Emission measure (f n e nu dV where d is the distance in cm and n e and nu 
are the electron and H densities in cm _3 ). c In units of total photons/cm 2 /sec in the line. d In units of photons/keV/cm 2 /sec at 1 
keV. e For the energy range 0.7-10.0 keV. The luminosity estimates are for an assumed distance of 2 kpc to HB3. 



